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摘 要 : 根据 收集 的 格尔木 地 区 地 下 水 位 动态 资料 .各 水 库 出 入 库 流量 资料 以 及 泉 集 河 实测 资料 等 ,基于 数值 模 
拟 方法 ,分 析 了 温泉 水 库 .南山 口 一 级 电站 以 及 洪水 等 因素 对 格尔木 河 地 下 水 位 上 升 程度 的 影响 。 结 果 表 明 : 格 尔 


木 市 地 下 水 位 的 上 升 是 南山 口 一 级 电站 洪水 以 及 温泉 水 库 调节 共同 作用 的 结果 。 南 山口 一 级 电站 的 蕾 水 使 得 市 


区 地 下 水 位 上 升 0.44 m; 洪 水 使 得 市 区 地 下 水 位 4.32 m; 温 凡 水 库 的 调节 作用 ,使 地 下 水 位 下 降 0.27 m; 地 下 水 的 大 


规模 开采 有 助 于 抑制 地 下 水 位 的 上 升 。 
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格尔木 是 青海 的 重要 城市 之 一 ,地 处 柴 达 木 盆 
地 。 全 国 最 大 的 干 润 内 陆 盐 湖 一 一 依托 察 尔 汗 盐 
湖 位 于 其 中 ,是 目前 全 国 最 大 的 钾肥 工业 生产 基 
地 。 格 尔 木 地 区 气候 干旱 ,但 受 格尔木 河 等 源 自 昆 
仑 山河 流 的 补给 作用 ,区 域内 地 下 水 资源 较为 丰 
富 。 地 下 水 位 周期 性 上 升 现象 持续 出 现 ,1977 年 以 
来 ,格尔木 河上 流 陆 续 修建 了 温泉 水 库 、 大 干 沟 水 
库 电 站 小 干 沟 水 库 电 站 、 乃 吉 里 水 库 电站 、 南 山口 
一 级 水 库 电站 等 水 利 设施 "。 随 着 水 利 工程 的 不 断 
建设 完善 以 及 大 洪水 的 发 生 , 近 20 a 来 格尔木 地 区 
地 下 水 位 处 于 历史 较 高 水 平 。2010 年 ,格尔木 市 出 
现 了 地 下 水 位 上 升 造成 的 严重 地 质 灾害 ,大 部 分 的 
农田 无 法 耕种 ,沼泽 湿地 面积 不 断 扩大 ,地 下 水 大 
量 泄 出 地 表 , 渍 水 泛滥 成 灾 ,致使 建筑 物 损 失 和 倒 
塌 , 工 厂 .学校 .企业 单位 .居民 点 .商场 等 面临 搬迁 ， 
直接 损失 达 6.0 x 10 元 左右 ,间接 损失 难以 估计 2 。 

长 期 的 地 下 水 位 抬升 不 仅 造成 城市 水 文 环 境 
的 变化 ,还 会 产生 其 他 潜在 的 危害 。 如 造成 土壤 的 
盐 渍 化 ,地 下 水 对 宕 土 及 建筑 物 的 腐蚀 增强 , 细 粉 
砂 液化 出 现 管 涌 等 现象 ,出 现 地 下 洞 室 基 础 上 浮 、 
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建筑 物 失 稳 等 现象 5”。 因 此 , 查 明 地 下 水 位 上 升 的 
原因 ,对 格尔木 地 区 的 地 质 灾害 预防 具有 至 关 重 要 
的 作用 。 已 有 的 研究 表明 :格尔木 地 区 的 地 下 水 位 
抬升 可 能 是 天 然 洪 水 与 水 利 工程 修建 共同 作用 的 
结果 2” 。 然 而 ,这 些 研究 只 提供 了 地 下 水 位 变化 与 
某 些 关键 因素 的 定性 关系 ,未 建立 机 理 模型 对 地 下 
水 位 致 灾 性 抬升 过 程 进行 定量 分 析 。 

本 文 的 研究 目标 ,是 建立 格尔木 地 区 的 地 下 水 
非 稳定 流 场 数值 模型 ,考虑 水 利 工程 的 作用 设置 不 
同 的 模拟 情景 ,评价 2009 一 2010 年 地 下 水 位 致 灾 性 
抬升 的 主要 成 因 , 用 于 预测 格尔木 地 下 水 位 的 变化 
趋势 ,并 用 2011 一 2018 年 地 下 水 位 的 相关 数据 进行 
了 验证 ,以 便 为 格尔木 市 预防 地 下 水 位 抬升 导致 的 
灾害 提供 决策 依据 。 


1 研究 区 概况 


格尔木 地 区 位 于 柴 达 木 盆地 南 缘 ( 图 1) , 深 居 
内 陆地 处 高 原 , 冬 季 寒 冷 ,夏季 凉 严 , 导 人 夜 温差 大 ， 
干燥 少雨 ,风速 强劲 ,沙尘暴 多 等 高 原 气 候 特 征 , 属 
典型 的 高 原 内 陆 盆 地 干旱 气候 。 主 要 河流 为 : 格 尔 
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图 例 E 铁路 | A Oe) 格尔木 市 政府 驻地 工作 区 范围 


图 1 研究 区 位 置 示 意图 


Fig. 1 Study areas location map 


木 东 河 、 格 尔 木 西河 ` 红 旗 河 、 金 水 河 ` 巴 水 河清 


数据 ;(2) 格尔木 市 水 库 管 理 处 提供 的 1980 一 2017 


水 河 。 年 力 吉 里 水 库 等 出 库 入 库 流量 资料 ; (3) 1992 一 
格尔木 地 区 山 前 冲 洪 积 平原 地 带 ,沉积 了 巨 厚 ” 1997 年 崇 达 木 综 合 地 质 勘查 局 多 次 实测 的 格尔木 


的 第 四 系 松散 沉积 物 ,为 地 下 水 赋 存 提供 了 良好 的 
空间 。 中 、 上 更 新 统 为 冰 磺 、 冰 水 、 冲 湛 积 地 层 , 岩 


市 渠 系 入 渗 量 ;(4) 格尔木 市 水 利 部 门 发 布 的 
2008 一 2017 年 渠 系 灌溉 量 ; (5) 2008—2017 年 格 尔 


性 主要 为 砂 卵 砾石 ` 泥 质 砂 砾石 ,岩层 颗粒 善 般 较 
粗 ,松散 ,相应 的 孔隙 度 、 给 水 度 也 大 ,在 接受 地 表 
水 的 垂直 渗 漏 补给 后 ,成 为 地 下 水 分 布 和 赋 存 的 主 
要 空间 。 地 下 水 类 型 主要 为 扇 区 的 第 四 系 松 散 兰 
类 孔 际 潜水 和 溢出 带 前 缘 的 少量 承 压 水 。 


木 市 主要 果 集 河 的 实测 资料 ;(6) 1990 年 以 来 , 格 尔 
木 市 地 下 水 位 长 观 孔 的 实测 数据 。 


3 ”原理 与 方法 
3.1 地 下 水 流 数值 模拟 方法 


2 数据 来 源 选用 瑞士 联邦 苏黎世 理工 大 学 开发 的 Process- 
ing Modflow Pro OEREN 了 区 域 地 下 水 非 稳定 流 
数学 模型 涉及 的 数据 主要 有 :(1) 青海 省 格 尔 ”的 数值 模拟 *%。 模 拟 范围 是 格尔木 河流 域 平原 区 ， 


木 市 水 文 站 提供 的 1959 一 2017 年 的 断面 流量 观测 


青海 省 环境 地 质 勘 查 局 


包括 地 下 水 流 场 模型 和 河流 - 果 集 河 模 型 和 %%。 以 


.青海 省 格尔木 河 冲 洪 积 遍地 下 水 资源 勘查 [R]. 2018. 


© 青海 省 水 文 地 质 工 程 地 质 勘 察 院 .青海 省 柴 达 木 盆地 地 下 水 资源 调查 [R]. 2002. 
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高 斯 投影 地 图 为 底 图 ,用 正方 形 网 格 对 建 模 区 进行 
剖 分 ,坐标 系 与 公里 网 平行 , 剖 分 为 1 kmxl km 正方 
形 网 格 ,与 地 图 公里 网 一 致 ,南北 向 前 分 101 格 , 东 
西向 62 格 ,有 效 面 积 $197 km。 分 别 用 平均 水 文 过 
程 和 典型 年 组 合 水 文 过 程 模 拟 4 种 方案 下 的 地 下 水 
变化 趋势 ,评价 2009 一 2010 年 格尔木 市 地 下 水 位 致 
灾 性 抬升 的 主要 成 因 。 模 型 的 应 力 期 (stress peri- 
od) 时 间 长 度 为 月 ,计算 的 时 间 步 长 为 天 。 
3.2 地 下 水 位 抬升 的 模拟 情景 方案 

为 分 析 2009 一 2010 年 地 下 水 位 上 升 机 理 , 弄 清 
各 因素 对 地 下 水 环境 影响 程度 , 设 定 5 种 情景 开展 
模拟 分 析 。 各 情景 方案 如 下 : 

方案 一 :该 方案 为 参照 背景 ,情景 中 无 温泉 水 
库 调 节 南山 口 一 级 电站 未 蕾 水 并 且 考 虑 乃 吉 里 电 
站 渗 漏 ,模拟 2008 之 前 地 下 水 状态 ,作为 分 析 其 他 

影响 因素 参照 背景 。 

方案 二 :专门 考虑 温泉 水 库 调节 对 地 下 水 的 影 
响 ,与 方案 一 相 比 ,将 格尔木 水 文 3 站 径流 过 程 替换 
为 温泉 水 库 调节 后 的 多 年 (1959 一 2007 年 ) 平 均 过 
程 。 其 他 模型 设置 与 方案 一 完全 相同 。 与 参照 背 
景 模拟 结果 相 比 较 , 可 分 离 温 果 水 库 调节 对 地 下 水 
环境 影响 。 

方案 三 :专门 考虑 南山 口 一 级 电站 水 库 渗 漏 
对 地 下 水 影响 ,在 方案 二 基础 上 ,增加 南山 口 一 级 
电站 渗 漏 ,其 他 与 方案 二 设置 相同 。 该 方案 与 参 
照 背景 方案 比较 ,可 分 离 温泉 水 库 调 节 、 一 级 电站 
渗 漏 共同 作用 对 地 下 水 的 影响 。 与 方案 二 模拟 结 
果 比 较 , 可 分 离 一 级 电站 渗 漏 对 地 下 水 环境 的 
影响 。 

方案 四 :考虑 2008 一 2012 年 连 丰 特大 洪水 影 
响 ,在 方案 三 基础 上 , 即 在 温泉 水 库 调 节 J3 EK 
库 与 南山 口 电站 水 库 渗 漏 基 础 上 ,接续 模拟 2008 一 
2012 年 河流 连 丰 \ 特 丰 环 境地 下 水 过 程 。 格 尔 木 水 
文 3 站 径流 过 程 取 2008 一 2012 年 3 站 流量 ,各 水 源 
群 开采 量 取 2008 一 2012 年 实际 开采 量 ,其 他 与 方案 
三 相同 。 该 方案 与 方案 三 结果 比较 ,可 分 离 连 丰 、 
特大 洪水 对 地 下 水 环境 的 影响 。 

方案 五 :将 2011 一 2018 年 实际 地 表 水 .水 库 实 
际 调 水 .地 下 水 开采 量 以 及 规划 水 源 情况 等 带 入 模 
型 ,进一步 验证 原 有 模型 并 且 模 拟 洪水 退却 时 水 位 
的 变化 情况 。 


4 各 模拟 方案 结果 对 比分 析 


4.1 方案 一 模拟 结果 

该 方案 模拟 结果 作为 其 他 方案 参照 背景 ,以 便 
分 离 各 种 因素 对 环境 的 影响 。 模 拟 水 均衡 分 析 结 
果 见 表 1, 浅 层 地 下 水 模拟 流 场 见 图 2 ,潜水 地 下 水 
位 埋 深 见 图 3。 

在 参照 条 件 下 ,格尔木 河 3 站 入 境 量 7.89x108 
ma (216.2x10 md), 其 中 ,河水 渗 漏 补给 5.49x 
10° m*+a'(150.3410! mq ) ,灌溉 引水 1.77x10° m+ 
a (48.1x10' md ), 两 者 消耗 占 入 境 流量 的 92% , Fl 
R 0.63x10° ma'(17.76x10: m- d’), HIERE 
河 溢出 量 1.02x10* m+ a!(27.95x10' m°: d), RIA 
1.65x10° m°- a!(45.15x10* mq") 经 格尔木 西河 流向 
东 达 布 逊 湖 ,其 中 部 分 消耗 于 沿途 蒸发 , 余 者 入 东 

在 格尔木 河 天 然 水 文 条 件 下 ,平原 区 多 年 平均 
地 下 水 总 补给 量 为 6.38x10' ma (174.88x10' md )， 
其 中 水 文 3 站 至 万 吉 里 水 库 段 渗 漏 补给 0.73x108 
ma (20.03x10' md 站 ), 占 总 量 的 11.5% , 乃 吉 里 水 
库 至 市 区 河 段 渗 漏 4.76x10: m-a”, 4 74.6% ,两 者 
之 和 为 河水 向 地 下 总 补给 量 5.49x10* m- a '(150.34x 
10 m°- d`’), Fy 86.1% , 即 地 下 水 补给 量 绝 大 部 分 由 
河水 转化 形成 ,其 他 补给 0.89x1l0xs m +a’, AX i 
13.9%. 

BRAC AWARE ATE RIK AZ RA 
与 补给 量 相对 应 的 排泄 特征 ,总 排泄 与 总 补给 量 相 
等 , 亦 为 6.38x10: ma.arl(174.88x10: m°- d7), HEPR 
(HM KHEN 2.33 10% m+ a°'(63.78x 10" m°: d’), 
36.5% , Z$ KE 3.24108 m+ a (88.74 10° m+ d`’), 
50.7% , Bil HF ZK HENE 2 KERI RE F RK i h AN 
RR HE FKF IR GE 0.37x 108 m-a” (10.25x10 m- d°), 
仅 占 5.8% ,向 北部 荒漠 (沙漠 . 盐 漠 ) 地 下 径流 0.44x 
108 ma.a(12.11x10:m.d)( 大 部 分 耗 散 于 北部 荒漠 
区 浅 埋 带 蒸发 ), 仅 占 6.93%。 

4.2 方案 二 模拟 结果 

本 方案 分 析 结 果 与 背景 方案 相 比 较 , 分 离 出 温 
泉水 库 调节 对 地 下 水 的 影响 。 

与 参照 背景 相 比较 ,经 温泉 水 库 调节 使 得 格 尔 
木 河流 量 在 枯 期 相对 增 大 ,汛期 相对 减少 ,改变 了 
年 内 径流 过 程 ,使 得 地 下 水 入 渗 量 在 年 内 朝 着 均一 


202012.00051v1 


chinaXiv 


汪 生 斌 等 :基于 数值 模拟 的 格尔木 地 区 地 下 水 位 致 灾 性 抬升 机 到 


时 


表 1 各 方案 均衡 要 素 


研 


ChinaXiv 合 作 期 刊 


Tab.1 Results of water balance for different scenarios /(10* m*+d") 
地 下 水 均衡 要 素 方案 一 方案 二 方案 三 方案 四 方案 五 
地 下 水 补给 格尔木 河沿 途 渗 漏 补给 130.51 128.75 111.13 169.49 133.09 
乃 吉 里 电站 水 库 渗 漏 20.03 20.03 20.036 20.03 20.29 
南山 口 一 级 电站 水 库 渗 漏 0 0 26.96 29.96 21.00 
肩 区 沟谷 渗 漏 补 给 9.46 9.46 9.46 9.46 10.77 
渠 系 灌溉 渗 漏 补给 14.78 14.78 14.78 14.78 16.79 
昆仑 山地 下 径流 补给 0.104 0.10 0.10 0.10 0 
补给 项 合计 174.88 173.12 182.46 243.82 202.08 
地 下 水 排泄 开采 地 下 水 10.25 10.25 10.25 14.3 100 
格尔木 河 与 录 集 河 溢出 63.78 63.54 65.59 77.55 59.86 
浅 埋 带 与 沼泽 湿地 蒸发 88.73 87.22 94.49 113.97 40.27 
向 北部 荒漠 区 径流 12.12 12.12 12.14 12.17 1.97 
排泄 项 合计 174.88 173.12 182.46 217.99 202.08 
地 下 水 补给 -排泄 均衡 差 0 0 0 25.83 0 
地 表 水 分 析 格尔木 河 (三 站 ) 入 境 量 216.23 216.23 216.23 382.08 223.7 
灌溉 引水 48.49 48.49 48.49 48.49 48.49 
格尔木 河 干流 沿 河 渗 漏 + 库 区 渗 150.53 148.77 158.11 219.48 174.38 
格 西河 及 汇 入 支流 溢出 63.78 63.54 65.59 77.55 76.32 
向 北部 荒漠 ( 盐 漠 ) 径 流 80.98 82.49 75.21 191.67 59.9 
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图 2 格尔木 河 位 置 示意 图 及 方案 一 模拟 的 潜水 等 水 位 线 


Fig. 2 Schematic diagram of the Golmud River and modeling results of phreatic water table for the scenario—1 
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图 3 溢出 带 附近 洪水 直 


WR TR) 


Fig. 3 Depth to water table in the near overflow zone (scenario~ 1) 


化 的 方向 发 展 ,但 未 改变 补给 总 量 。 水 库 调节 后 ， 
河流 总 渗 汤 量变 化 不 多 ,上 略微 减少 , 果 水 和 地 下 水 
蒸发 亦 咯 有 下 降 ( 表 1)。 尽 管 总 渗 漏 量变 化 不 多 ， 
由 于 村 期 流量 增 大 ,流量 随 季 市 变化 更 加 平稳 , 治 
河 地 带 地 下 水 位 有 一 定 变化 ,市 区 地 下 水 位 平均 下 
降 0.27 m, 年 内 水 位 变 幅 略 减 小 ,由 原 2.07 m 降 为 
1.75 m, 由 此 可 见 , 温 果 水 库 的 调节 作用 使 得 地 下 水 
水 位 不 升 反 降 , 故 而 温泉 水 库 调节 不 是 诱发 市 区 水 
位 上 升 因 素 。 

4.3 方案 三 模拟 结果 

本 方案 计算 结果 与 背景 方案 .水库 调 节 方 案 比 
较 , 可 分 离 出 一 级 电站 渗 漏 对 地 下 水 环境 影响 。 

本 方案 水 均衡 分 析 列 于 表 1。 一 级 电站 水 库 渗 
漏 影响 地 下 水 位 变 幅 等 值 线 见 图 4。 

据 模 拟 分 析 ,2008 年 南山 口 一 级 电站 蓄 水 后 ， 
电站 水 库 大 量 渗 漏 ,河流 总 渗 漏 量 增加 ,与 方案 二 
相 比 , 渗 漏 量 与 渗 漏 位 置 都 发 生 了 明显 改变 。 格 尔 
木 河 上 段 增加 电站 渗 漏 26.96 x 10'm?+d"(3.12 m-s”), 


图 4 南山 口 一 级 电站 诱发 的 地 下 水 变 幅 等 值 线 
Fig. 4 Increment of groundwater level induced by the first 


Nanshankou power station 
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导致 下 游 河 水 流量 有 所 降低 ,使 格尔木 河 下 段 渗 漏 
减少 17.62x10' md ', 抵 前 后 河流 渗 漏 总 量 增加 
9.34x10' mx…d ,诱发 地 下 水 流 场 .市 区 地 下 水 位 、 泉 
水 流量 等 地 下 水 环境 发 生变 化 。 

在 格尔木 河 洪 积 扇 调节 作用 下 ,一 级 电站 渗 漏 
对 环境 影响 具有 滞后 性 ,形成 稳定 影响 需 清 后 10 a 
时 间 ( 图 5、 图 6)。 以 市 区 水 位 变化 为 例 , 一 级 电站 


蓄 水 发 电 前 相 比 ,影响 至 稳定 后 ,市 区 水 位 升 高 
0.44 m( 图 4)。 

南山 口 一 级 电站 水 库 蓄 水 后 ,对 水 环境 影响 具 
有 了 明显 滞后 性 ,滞后 期 3 a, 达 稳定 需 10 a; 电站 水 库 
B iii 26.96 10° md ,使 格尔木 河 下 段 渗 漏 减 少 
17.62x10' md ,两 者 抵 削 后 格尔木 河 渗 漏 量 净 增 
9.34x10' m’…d-'; 泉 集 河 溢出 增 量 2.06x10: md’ NT 
区 地 下 水 位 比划 水 发 电 前 平均 升 高 0.44 m, 对 地 下 


蕾 水 初期 , 因 河 流下 段 渗 漏 量 减少 ,市 区 水 位 不 仅 
不 升 ,反而 下 降 ,3 a 后 恢复 到 鞭 水 前 地 下 水 位 ,经 组 
慢 调 整 ,5 a 后 上 升幅 度 达 到 稳定 影响 值 的 63% ,10 a 
后 水 位 影响 达到 稳定 。 泉 集 河 灌 后 溢出 特征 与 水 
位 过 程 相 似 。 电 站 水 库 2008 年 8 月 开始 蓄 水 , 按 3a 
滞后 期 估算 ,至 2011 年 (3 a) ,水 库 渗 漏 影响 才刚 刚 
显现 。 


15 20 25 30 
a 0.20 时 间 /a 


图 5 市 区 地 下 水 位 对 南山 口 一 级 电站 的 滞后 响应 过 程 


Fig.5 The delayed response of groundwater level in the urban 


area to the first Nanshankou power station 
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图 6 南山 口 一 级 电站 水 库 渗 漏 , 泉 集 河 溢出 量 
潍 后 调整 过 程 


Fig.6 Nanshankou first-level power station reservoir leakage, 


Quanji river overflow lag adjustment process curve 


电站 水 库 渗 漏 对 区 域 地 下 水 流 场 影响 :水 库 渗 
漏 影响 至 稳定 后 ,与 方案 二 相 比 ,电站 水 库 附 近 形 
成 水 位 上 升 中 心 , 上 升幅 度 约 20 m, 上 升 范围 波及 
整个 肩 区 ,北部 延伸 至 泉水 溢出 带 ( 图 5)。 影 响 稳 
定 后 , 果 集 河 排 汇流 量 增 大 , 增 量 2.05x10'm…d ;与 


水 环境 有 一 定 影响 。 
4.4 方案 四 模拟 结果 

本 方案 以 2008 年 流 场 作为 初始 条 件 ,格尔木 河 
流量 过 程 修改 为 2008 一 2012 年 实际 连 丰 水 文 过 程 ， 
一 级 电站 水 库 2008 年 8 月 蓄 水 ,各 水 源 群 取 2008 一 
2012 年 实际 开采 量 。 

2008 年 格尔木 河 径 流量 9.05x10 ma 
(247.95x10* m°: d`) , 属 偏 丰年 (P=25% ,已 指 河水 流 
量 ) , 2009—2012 年 ,格尔木 河 径流 量 为 10.95x10*、 
18.95x 10°, 12.45x10°,13.43x10° ma ,连续 4 a 都 超 
过 二 十 年 一 遇 的 丰 水 年 ,其 中 ,2010 年 为 千年 一 遇 
特大 丰 水 年 。2009 一 2012 年 水 均衡 分 析 、 与 各 方案 
对 比 结果 见 表 1。 

对 比 不 同方 案 果 集 河 溢出 增 量 :2009 一 2012 
年 ,在 考虑 清 后 因素 情况 下 ,一 级 电站 蓄 水 (方案 
三 ) 泉 集 河 溢出 水 量 , 比 蓄 水 前 (方案 二 ) 仅 增加 0.7x 
10° md ', 而 特 丰 、 连 丰年 泉 集 河 溢出 水 量 比 方案 
二 增加 13.32x10' m°- d> ,比方 案 三 增加 12.62x10': m°- 
d, BRE 、 连 丰年 录 集 河 增 量 超过 一 级 电站 蕾 水 
影响 6 倍 以 上 。 各 主要 因素 对 地 下 水 位 抬升 至 灾 的 
影响 程度 :在 考虑 滞后 因素 情况 下 ,至 2012 年 ,一 级 
电站 蓄 水 引起 市 区 水 位 升幅 0.44 m ,而 洪水 引起 水 
位 升幅 4.32 m, 即 洪水 发 生年 对 市 区 水 位 升幅 的 影 
响 超 过 一 级 电站 渗 漏 8 倍 以 上 。 

通过 对 比 泉 集 河 溢出 量 、 市 区 水 位 升幅 表明 ， 
2010 一 2012 年 市 区 水 位 上 升 、 果 集 河流 量 增加 等 地 
下 水 环境 变化 ,主要 是 格尔木 河 特 丰 、 连 丰 造 成 的 ， 
昌 一 级 电站 鞭 水 有 一 定 的 支持 作用 ,与 格尔木 河 特 
丰 、 连 丰 相 比 很 弱 。 

为 考察 2012 年 以 后 地 下 水 演化 趋势 ,本 方案 预 
测 时 间 延 长 5a,2012 年 后 河水 流量 , 设 定 为 多 年 平 
均值 。 采 用 数学 模型 对 地 下 水 位 立 值 进行 模拟 "*， 
结果 表明 ,假设 后 续 为 平水 年 ,3 a 后 ,市 区 地 下 水 位 
15 25 Tt Te Ye E (EL 7K iz pv Fi 2805.5 m 以 下 (图 7)， 
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2010 年 洪水 果 。 南 山口 一 级 电站 革 水 使 得 泉 集 河 流量 增 大 0.7x 
2807 
10° md! ,市 区 地 下 水 位 上 升 0.44 m; 连 续 特 丰年 产 
生 的 洪水 使 得 泉 集 河流 量 增 大 13.32x10' m°- d, TH 
临界 水 位 值 区 地 下 水 位 4.32 m; 温 泉水 库 起 调节 作用 ,使 地 下 水 
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图 7 方案 四 模拟 的 市 区 地 下 水 位 长 期 变化 


Fig.7 Modeling results of the long-term groundwater 


1000 


level variation in the urban 


预测 的 整体 变化 趋势 与 图 4 实测 数据 基本 一 致 。 
4.5 方案 五 模拟 结果 

本 方案 在 方案 四 的 基础 上 将 时 间 序 列 延 长 至 
2017 年 ,并 将 水 源 开 采 量 调整 为 已 批复 的 格尔木 河 
冲 洪 积 遍地 下 水 开采 量 100x10* m°. d'o 

该 方案 下 开采 地 下 水 100x10:m.d ,地 下 水 环 
境 有 明显 变化 ,水 位 埋 深 增 大 ,泉水 .蒸发 减少 ,向 
下 游 盐 湖区 水 流 减 少 。 开 采 水 量 95% 来 自 夺 取 泉 
水 和 蒸发 ,泉水 由 4.4x10* m°- a™(120.6x10* m°- do") ak 
> 28 2.19108 ma '(59.9x10’ m°: d); 48 KH 2.23x 
10° m°- a (61.2x10 m- TAW 1.47 10° ma 
(40.3x10' m°- d'o RIK ite AAS A a REE 
接受 范围 。 开 采 后 格尔木 西河 和 果 集 河 减少 入 湖 
水 量 1.8x107 ma 。 对 盐湖 矿区 采 讽 等 用 水 有 一 定 
影响 ,建议 通过 那 棱 格 勒 河 调 水 工程 (已 竣工 ) 
缓解 。 

综 上 所 述 ,格尔木 河 特 丰 、 连 丰年 洪水 ,是 引起 
市 区 及 其 北部 地 下 水 位 上 升 的 决定 因素 ,南山 口 一 
级 电站 水 库 滩 漏 起 次 要 作用 。 温 果 水 库 调节 对 格 
尔 木 市 区 的 地 下 水 位 抬升 没有 影响 。 


5 结论 


针对 格尔木 市 地 下 水 位 上 升 至 灾 的 原因 ,本 文 
从 地 下 水 数学 模型 的 角度 ,通过 对 可 能 的 因素 进行 
了 分 析 , 得 出 以 下 结论 : 

(1) 格尔木 市 地 下 水 位 的 上 升 是 南山 口 一 级 电 
站 、2010 年 大 洪水 以 及 温 果 水 库 调节 共同 作用 的 结 


位 向 稳定 方向 发 展 ,该 水 库 修 建 后 市 区 地 下 水 位 下 
降 0.27 m; 地 下 水 开采 量 达 到 规划 的 100x10' m -d~ 
时 , 泉 集 河 洒 出 量 减 少 了 60.7x10¢ m -d ,有 助 于 抑 
制 地 下 水 位 的 上 升 。 

(2) 格尔木 地 区 地 下 水 位 上 升 机 理 : 由 于 2008 
年 以 来 格尔木 河 连续 丰 水 年 ,特别 是 2010 年 发 生 的 
大 洪水 和 水 库 渗 漏 作用 ,改变 了 格尔木 河 在 冲 洪 积 
扇 上 的 流 态 ,增强 了 河流 主要 是 下 河 段 河 水 渗 漏 强 
度 , 增 大 了 地 下 水 补给 的 增加 量 ,因而 也 增 大 了 地 
下 水 储存 量 的 增加 量 , 有 限 空 间 内 储存 量 的 增加 必 
然 会 导致 地 区 地 下 水 位 剧烈 上 升 。 

(3) 利用 数学 模型 预测 了 5 种 不 同 的 方案 , 认 
为 本 次 地 下 水 位 上 升 的 主要 原因 是 2008 年 以 来 格 
尔 木 河 连 续 丰 水 年 ,特别 是 2010 年 发 生 的 大 洪水 造 
成 的 ,南山 口 一 级 水 库 渗 漏 所 起 作用 比 前 者 弱 很 
多 ,大 规模 开采 地 下 水 有 助 于 抑制 地 下 水 位 的 上 
升 ,而 温泉 水 库 对 地 下 水 位 上 升 基本 无 影响 。 

(4) 数学 模型 预测 结果 表明 : 按 2012 年 后 按照 
平水 年 设 定 输入 条 件 ,3 a 后 地 下 水 位 将 逐渐 回落 至 
临界 水 位 以 下 ,可 缓解 水 位 上 升 灾情 。 
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Study on mechanisms of geological hazards caused by groundwater level 


rising in the Golmud area based on numerical simulation 
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Abstract: From 2009 to 2010, the groundwater level in Golmud continued to rise significantly, the area of marshes 


and wetlands continued to expand, and a large amount of groundwater leaked out of the surface, resulting in 


waterlogging and flooding. As a result, buildings were damaged and collapsed, farmland was flooded, shopping 


malls were flooded, and some bungalows caused cracks or collapses in the walls. Residential buildings were greatly 


threatened. All of these caused great economic losses. In this paper, based on the collected long series of 


groundwater level dynamic data, the inflow and outflow data of various reservoirs and the measured data of the 


Quanji River combined with a mathematical model were used to analyze the causes of the rise of groundwater level. 


It was found that the direct cause of the groundwater level rise that led to disaster was flooding. 


Keywords: groundwater; level rise; numerical simulation; flood; Golmud area 


